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Quantum error correction is required to compensate for 
the fragility of the state of a quantum computer. We report 
the first experimental implementations of quantum error cor- 
rection and confirm the expected state stabilization. In NMR 
computing, however, a net improvement in the signal-to-noise 
would require very high polarization. The experiment imple- 
mented the 3-bit code for phase errors in liquid state state 
NMR. 

PACS numbers: 03.65.Bz, 89.70.+c,89.80.th,02.70.-c 



Quantum computers exploit the superposition prin- 
ciple to solve some problems much more efficiently 
than any known algorithm for their classical counter- 
parts. These problems include factoring large numbers 
||J (thereby breaking public key cryptography), combina- 
torial searching Q and simulations of quantum systems 
|^-^| . Exploiting the power of quantum computation was 
thought to be physically impossible due to the extreme 
fragility of quantum information This judgment 

was demonstrated to be overly pessimistic when quantum 
error-correction techniques p|-|l(]|] were found to protect 
quantum information against corruption due to imper- 
fect control and decoherence. It is now known that for 
physically reasonable models of decoherence a quantum 
computation can be as long as desired with arbitrarily ac- 
curate answers, provided the error rate is below a thresh- 
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Thus decoherence and imprecision are 
: ered insurmountable obstacles to realiz- 
ing a quantum computer. 

The chief remaining obstacle to quantum computing is 
the difficulty of finding suitable physical systems whose 
quantum states can be accurately controlled. Devices 
based on ion traps ]ll| have so far been limited to two 
bits [|l6j. Recently, liquid state NMR techniques have 
been shown to be capable of quantum computations with 
at least three bits |l7],[R|. This makes it possible, for 
the first time, to experimentally implement the simplest 
quantum error-correcting codes, and so test these ideas in 
physical systems with a variety of decoherence, dephasing 
and relaxation phenomena. 

In room temperature liquid state NMR, one can coher- 
ently manipulate the internal states of the coupled spin 
^ nuclei in each of an ensemble of molecules subject to 
a large external magnetic field. Although the ensemble 
nature of the system and the small energy of each spin 
imply that the set of accessible states is highly mixed, it 
has been shown that experimental methods exist that can 



be used to isolate the pure state behavior of the system, 
thus permitting limited application of NMR to quantum 
computation p9| , p0[ . 

Here we describe the implementation of a quantum 
error-correcting code which compensates for small phase 
errors due to fluctuations in the local magnetic field. The 
behavior of this code was measured for two systems: The 
13 C labeled carbons in alanine subject to the correlated 
phase errors induced by diffusion in a pulsed magnetic 
field gradient, and the proton and two labeled carbons 
in trichloroethylenc (TCE) subject to its natural relax- 
ation processes. In alanine, we observed correction of 
first-order errors for a given input state, while in TCE 
we demonstrated state preservation of an arbitrary state 
by error correction. 

In NMR loss of phase coherence may be both due to 
bona fide irreversible decoherence |2l| and due to de- 
phasing which can be reversed by spin echo. Both lead 
to symptoms which, in an ensemble, are identical from 
the standpoint of quantum computation. Both can be 
corrected using the same scheme we shall demonstrate 
below. 

It is important to comment on the applicability of this 
technique to ensemble quantum computing. Although 
our experiments validate the usefulness of error correc- 
tion for quantum computing with pure states, there is a 
substantial loss of signal associated with the use of an- 
cilla spins in weakly polarized systems. We argue that 
in this setting, the loss of signal involved in exploiting 
ancillas removes any advantage for computation gained 
by error correction, at least unless the system is suffi- 
ciently polarized to enable the generation of nearly pure 
states. Nevertheless, our experiments demonstrate that 
error-correcting codes can be implemented, and that they 
behave as predicted. 

The simple three-bit quantum error-correcting code 
used here is designed to compensate to first order for 
small random phase fluctuations. These fluctuations con- 
stitute a random evolution of the state 
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where bi is or 1, Oi is a random phase variable, and <r l z is 
the Pauli matrix acting on the i'th spin. The Oi depend 
on the error rates in the model, which is described in 
detail below. 

The error-correcting code itself is a phase variant of 
the classical three bit majority code with a decoding 
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technique that preserves the quantum information in the 
encoded state. Let |+) = (|0) + \l))/\/2 and |-) = 
(|0> - |l))/\/2- The state (a|000) + 01 100)) is encoded 

as a\ + + +) + f3\ — ) by a unitary transformation. 

The first-order expansion of the operator in Eq. ([!]) in 
the small random phases is 
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which evolves the encoded state to 
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field. This is achieved by turning on an external field 
gradient V Z B = dB z /dz across the sample for a time 
5. This modifies the magnetization in the sample with 
a phase varying linearly along the z direction accord- 
ing to d<f>/dz — nS'ydBz/dz, where n is the coherence 
order of the density matrix element and 7 is the gyro- 
magnetic ratio. A reverse gradient is used to rcfocus 
the magnetization after allowing molecular diffusion to 
take place for amount of time t. As a result of random 
spin displacement Az, the phases of the spins are not re- 
turned to their original values but are randomly modified 
by (n5^dB z /dz)Az. For a gaussian displacement profile 
with a width of \/2Dt, the effective decoherence time of 
this process is proportional to the diffusion constant D 
as well as to the square of the coherence order n p3] : 



The different errors map the encoded state into or- 
thogonal subspaces. Thus one can determine which error 
occurred without destroying the encoded quantum infor- 
mation. This is accomplished by a measurement which 
reveals the subspace the state has moved into. After 
decoding, the original state of the first spin can then be 
restored by a unitary transformation, while the other two 
spins contain information (the syndrome) about the er- 
ror which occurred. A network which accomplishes the 
encoding, decoding and error-correction steps is shown 
in Figure |l|. 

In NMR experiments, non-unitary processes are classi- 
fied as spin-lattice and spin-spin relaxation p2| , p3| . The 
former involves an exchange of energy with the non-spin 
degrees of freedom (the lattice) , and returns the spins to 
thermal equilibrium. The latter leads to a loss of phase 
coherence. For spin ~ nuclei, both processes are due to 
fluctuating local magnetic fields. The three spin code 
corrects for errors due to locally fluctuating fields along 
the z axis. 

We focus on a weakly coupled three-spin system where 
the strongest contribution to coherence loss is from ex- 
ternal fields which contribute the Hamiltonian 

-H R = j 1 ! 1 ■ B\t) + 7 2 I 2 • B 2 (t) + 7 2 I 3 • B 3 (i) , (4) 

where I = (I x ,I y ,I z ) and I u = \<J U (u = x, y,z). In the 
case of slowly varying external fields the induced ran- 
dom phase fluctuations are identical to those described 
in Eq. (jl]). As a result, the off-diagonal elements of the 
density matrix decay exponentially at a rate which de- 
pends on the fields B fc at each spin, their gyromagnetic 
ratios 7 fc , the coherence order and the zero frequency 
components of the spectral densities of the fields. The 
"coherence order" is the difference between the total an- 
gular momenta along the z-axis of the two states \b), \b') 
(in units ofh/2) @. 

To obtain a clean demonstration of error correction, 
a simple error model was implemented precisely in the 
case of alanine. This implementation used the random 
molecular motion induced by diffusion in a constant field 
gradient to mimic the effect of a slowly varying random 
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This artificially induced "decoherence" in the alanine ex- 
periments is an example of completely correlated phase 
scrambling. This occurs naturally if all the spins have 
equal gyromagnetic ratios in the slow motion regime. 
We used TCE to demonstrate error-correction in the 
presence of the natural decoherence and dephasing of a 
molecule where T2 <C T\ . 

Most NMR experiments are described using the prod- 
uct operator formalism 



25Q. This formalism describes 
the state as a sum of products of the operators ij; , I* , ij. 
The identity component of such a sum is the same for any 
state and is usually suppressed to yield the "deviation" 
(traceless) density matrix. The effect of error-correction 
can be understood from the point of view of this for- 
malism. As an example, consider encoding the state \\ 
using two ancillas initially in their ground states. Up to 
a constant factor, the initial state is described by 
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The density matrix p A consists of an incoherent sum of 
four terms. 

After encoding the state is 



p B = + Il + I 3 x + 41^13) 
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In the case of completely correlated phase errors, this 
decays as 
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Decoding and error correction mixes these states together 
so as to cancel the initial decay of the first spin. The 
reduced density matrix for the first spin becomes 
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The effect of error correction can be seen from the ab- 
sence of terms depending linearly on t. Similar results 
can be obtained for other initial states and also for un- 
corrected phase errors. 

In the alanine experiments, each of the four product 
operators in the sum of Eq. | was realized in a sepa- 
rate experiment, and the final state after encoding and 
decoding inferred by adding the results. The loss of po- 
larization over time in each product operator was mea- 
sured explicitly in each experiment. The results were 
added computationally to simulate the effect of the Tof- 
foli gate and are shown in Figure |[ The components of 
I^I^Iz which evolved as a single and a triple coherence 
when encoded were separated by gradient labeling Q. 
The initial slopes of the decay curves for each opera- 
tor were estimated by linear square fit to the logarithm 
of the measured intensities. The values obtained are 
-0.0034 ( I*), -0.0030 (21*12), _ .0047 (21^), -0.0038 
(4I^lf, single coherence) and -0.0395 (4I^I|, triple 
coherence) When the slopes are added as required for 
error-correction with a pseudopure input, we get a mea- 
sured net slope of 0.00081, which should be compared 
to the slopes for the single coherences whose average is 
—0.0037. The net curve has quadratic behavior for small 
delays to within experimental error. 

The goal of our experiments with TCE was to establish 
the behavior of encoding/decoding and error correction 
on all possible initial states subject to the natural deco- 
herence and dephasing. The spins were prepared in the 
states 

P(il + I2)(il + I3), (10) 

with p one of the four inputs |l, 1\ 1 I*, I*. Any possible 
input is just a linear combination of these four states. 
We used gradient methods to directly generate the four 
states of Eq. [l^ with a series of experiments. They were 
then subjected to pulse sequences for encoding and de- 
coding (experiment I) with a variable delay between the 
two operations. Decoherence and dephasing take place 
during the delay. The reduced density matrix on the 
first spin (the output) was measured. In the second ex- 
periment (II) decoding was followed by error correction 
before the output was determined. Ideally the output 
would be identical to the input. The measured outputs 
were compared to the ideal ones by computing the "en- 
tanglement fidelity" ^f|. This is a useful measure of 
how well the quantum information in the input is pre- 
served. Entanglement fidelity is the sum of the correct 
polarization left in the output state for each input. More 
precisely, given input I*, let f a be the relative polariza- 
tion of I* in the output compared to the input. Then 
/ = |(1 + f x + f y + f z ), this formula is correct for pro- 
cesses which do not affect the completely mixed state \ 1 . 



The results for nine different delays are shown in Fig- 
ure ^|. The curves show that error correction decreases 
the initial slope by a factor of ~ 10 (by square fit to the 
logarithm). 

Our demonstration of error-correction does not imply 
that error-correction can be used to overcome the prob- 
lems of high temperature ensemble quantum computing. 
In this model of quantum computing, the initial state can 
be described as a small, linear deviation from the infinite 
temperature equilibrium. Thus, the deviation is propor- 
tional to a Hamiltonian of n weakly interacting particles. 
In this limit no method of error-correction based on ex- 
ternally applied, time dependent fields can improve the 
polarization of any particle by more than a factor propor- 
tional to ^/n [^7j . If one wishes to use error-correction an 
even bigger problem is encountered: The initial state of 
the ancillas used for each encoding/decoding cycle must 
be pure. In the high temperature regime, the best we can 
do to ensure that is to generate a pseudopure deviation 
in the ancillas. Unfortunately, this deviation has to be 
created simultaneously on all ancillas, leading to an expo- 
nential reduction in polarization proportional to the total 
number of ancillas required pi| . This always exceeds the 
loss of polarization. Another problem is the inability to 
reuse ancilla bits. This has two consequences. The first 
is that decoherence rapidly removes information in the 
state, leading to computations which are logarithmically 
bounded in time f29f| . Second, the total number of ancil- 
las required is proportional to the time-space product of 
the computation, rather than to a power of it's logarithm. 

Our work shows that liquid state NMR can be used to 
test fundamental ideas in quantum computing and com- 
munication involving non-trivial numbers of bits subject 
to a variety of errors. Our experiments demonstrate for 
the first time the state preserving effect of the three bit 
phase error-correcting code. The first-order behavior was 
established to high accuracy for a specific state in alanine, 
while the overall effect was observed and the improve- 
ment in state recovery verified in TCE. These experi- 
ments confirm not only the validity of theories of quan- 
tum error correction in a simple case, but also demon- 
strate the ability, in liquid state NMR, to control the 
state of three spin-half particles. This is an important 
advance for quantum computing, as this is the first sys- 
tem where this degree of control has been successfully 
implemented. 
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FIG. 1. Network for encoding, decoding and error correc- 
tion. The circuit describes the evolution of the 3 bits as a 
function of time. The gate «S*x corresponds to a control-not, 
i.e a gate which flips the target bit (x) if the control bit (•) 
is in the state [1). ^(90) represent a single bit gate corre- 
sponding to a rotation by an angle of n/2 around the y-axis. 
The Tofolli gate flips the target bit (x) if the two control bits 
(•) are in the state |1). A detailed implementation of these 
gates is given in [17]. The information carrying bit is carbon 
1 (see Figures 2 and 3) in both experiments. 
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FIG. 2. The intensities of the magnetization of the first 
spin after applying the dephasing and decoding procedures 
described in the text, together with single exponential fits 
to the intensities versus the dephasing time r. The three 
mixed states I*, I^I^, ijlz, evolved as single quantum co- 
herences during r, whereas I^IzIz evolved as a mixture of 
single and triple quantum coherences, which have been plot- 
ted separately (single and triple). The sum of these inten- 
sities and the corresponding fits (Error corrected) give the 
intensities and fit characteristics of the same experiment us- 
ing a pseudopure state (see text). The initial slope of the 
total is close to zero, thus showing that the error-correction 
procedure was able to cancel loss of coherence of the 
pseudopure state to first order. 
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FIG. 3. Experimentally determined entanglement fidelities 
for the TCE experiments after decoding (red) and after de- 
coding and error correction (green). The continuous curves 
are interpolations of the data points. The broken curves were 
determined by simulating the pulse sequence using the mea- 
sured coupling constants and estimated T2's of 1.1s (CI), 0.6s 
(C2) and 3s (H). Differences between experimental and the- 
oretical curves are attributed to lack of precise knowledge of 
the error model. Errors in the data points are approximately 
0.05. Note that since the proton T2 is much longer than that 
of the carbons, the long term gain in fidelity is partially due to 
recovery of polarization from the proton. The demonstration 
of error correction lies in the initial slope. The curves show 
that error correction decreases the initial slope by a factor of 
~ 10 (by square fit to the logarithm). 
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